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TURBULENT MIXING OF COAXIAL COMPRESSIBLE 
HYDROGEN-AIR JETS 

By James M. Eggers 
Langley Research Center 

SUMMARY 

An experimental and analytical study of the compressible turbulent mixing of par- 
allel coaxial hydrogen-air jets has been conducted. Data were acquired for outer air jet 
Mach numbers of 1.32 and 2.50. The inner hydrogen jet Mach number was approximately 
0.9 for both air jet Mach numbers. All jets had a total temperature near 300 K and mixed 
in an unconfined region at 1 atmosphere. Experimental hydrogen mass fraction profiles 
and velocity profiles were determined throughout the near-field and far -field mixing 
regions. 

The validity of different eddy viscosity models was studied by incorporating them 
into a finite-difference- type analysis and attempting to compute the present hydrogen-air 
data and previous air-air data. A formulation of eddy viscosity based on a mass flow 
defect (or excess) across the mixing zone was unsatisfactory in computing the data. The 
kinematic eddy viscosity model of Cohen and a kinematic form of an eddy viscosity model 
used in a previous study both satisfactorily correlated the data. It was concluded that a 
kinematic form of eddy viscosity which provides radial as well as axial variation in 
dynamic eddy viscosity, through incorporation of the local density, is essential in order 
to compute both hydrogen-air and air-air mixing data. 

INTRODUCTION 

An interest in the computation of turbulent mixing exists because of the wide variety 
of applications. Typical problems involving turbulent mixing occur in jet-engine exhaust- 
noise generation, shear-layer interference heating, ejector design, mixing of pollutants 
with the atmosphere, and fuel injector design. The latter problem of fuel injector design, 
specifically the mixing of hydrogen fuel with air in a supersonic-combustion ramjet- engine 
combustor, is the motivation for this study. The study is restricted to consideration of 
the unconfined mixing of nonreactive, circular coaxial jets, which is an approximation to 
the downstream parallel injection of fuel from an in-stream injector. 

In order to design a combustor for a hydrogen-fueled hypersonic ramjet engine, it 
is necessary to predict the fuel distribution obtained from a given injector design. For 



combustors of small height, it is conceivable to inject the fuel from the wall and achieve 
the desired uniform fuel distribution. As combustor size increases, however, it becomes 
necessary to consider injection from struts or, if the engine geometry permits, from an 
inlet center body in order to uniformly distribute the fuel. It is advantageous to inject the 
fuel in a downstream direction in order to benefit from the momentum of the fuel and in 
order to minimize pressure disturbances initiated by the fuel injection. Therefore, an 
analytical method is needed to predict the turbulent mixing of parallel compressible 
streams of hydrogen and air. 

Unfortunately, current knowledge of the fundamental nature of turbulence is not suf- 
ficient to permit generation of a completely analytical solution to the mixing problem. 
Analysis of the mixing of turbulent flow fields has been performed by employing the 
boundary-layer equations with suitable models for the turbulence terms (conventional eddy 
viscosity and turbulent Prandtl and Lewis numbers). The semiempirical nature of the 
analysis arises through the need to specify the magnitude of the turbulent Prandtl and 
Lewis numbers and an eddy viscosity model. No universally accepted eddy viscosity 
model is currently available. 

Many different eddy viscosity models have been proposed and reported in the litera- 
ture. In reference 1, Ferri proposed a model relating the dynamic eddy viscosity to a 
mass flux difference. Later studies reported in reference 2 noted no tendency for two jets 
to remain segregated (unmixed) when the mass flux difference approached zero as Ferri 's 
model would indicate. Also, no tendency was reported for two jets to remain unmixed 
when the velocity difference approached zero. From the studies of reference 2 it was not 
possible to determine whether the inadequacy of the Ferri model (or the inadequacy of 
Prandtl's model which expressed eddy viscosity as proportional to a velocity difference) 
was merely a limitation of the range of applicability of the models or an indication of 
invalidity. In reference 2, Alpinieri presented a model which was developed from purely 
empirical considerations and which correlated his data. 

A more recent study by Schetz (ref. 3) expressed the dynamic eddy viscosity as being 
proportional to the mass flow defect (or excess) across the mixing zone. Also given in 
reference 3 were the computations resulting from using the Schetz model, the Prandtl 
velocity difference model, and the Ferri model. The general observation reported in the 
reference was that the Schetz model provided results in better agreement with the data 
considered than any other previously suggested model. 

A kinematic eddy viscosity model has been developed by Cohen and reported in ref- 
erence 4. The Cohen model was satisfactorily applied to hydrogen-air, hydrogen-nitrogen, 
reacting, and nonreacting data (ref. 5). Comparison calculations using the Prandtl velocity 
difference model, the Ferri model, and the Cohen model for high-temperature hydrogen- 
nitrogen mixing were also presented in reference 5. It was noted that the Prandtl model 
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produced a much too slow mixing rate and the Ferri model produced a much too rapid 
mixing rate for the data considered. 

Previous to the present study the author had participated in a coaxial air-air mixing 
study (ref. 6). Initial efforts to correlate the data were made with an eddy viscosity model 
attributed by the authors of reference 7 to Zakkay. It was found that limitations imposed 
by the definition of the mixing width in the Zakkay model restricted its application to the 
far-fi^ld mixing zone. Therefore, a new mixing width was defined and resulted in 
improved data correlation and a new dynamic eddy viscosity model. 

] 

Some additional eddy viscosity models for the region downstream of the potential 
core and wake flow have been summarized in reference 3. Despite the sizable number of 
mixing studies performed, illustrated by the numerous suggested eddy viscosity models, 
these have not resulted in guidelines which could be expected to lead to the development 
of a satisfactory eddy viscosity model. Progress in developing a satisfactory model has 
been hampered by the limited quantity of data each investigator has considered as well as 
the insufficiencies and uncertainties of the available data. 

Data resulting from experiments involving nonreacting coaxial hydrogen-air mixing 
have been presented in references 1, 2, 5, 8, and 9. However, only references 1, 5, and 8 
contain data where the air stream was supersonic, which is of primary interest for 
supersonic- combustion- ramjet application. Of these three references, only reference 5 
(data referred to in ref. 5 are actually for hydrogen mixing with a high-temperature viti- 
ated nitrogen stream) gives detailed survey data for the initial mixing conditions. Even 
in the data of reference 5, initial boundary layers are not well defined because the data 
were taken at an elevated temperature. As a result, the size of the total temperature 
probe limited the resolution of the data. The supersonic data of references 1 and 8 dealt 
only with the far field, that is, the flow field downstream of the potential core. It has been 
noted in reference 10 that the near field has received less attention than the far field 
although the near field is considered to be just as important in the overall design of a 
ramjet combustor. In summary, available supersonic hydrogen-air mixing data do not 
adequately define initial conditions or the near -field mixing region. 

The purpose of the present investigation was twofold. The first objective was to 
generate detailed supersonic coaxial hydrogen-air mixing data, both in the near field and 
far field, which could be used to aid in the development of an analysis and eddy viscosity 
model. The second objective was to evaluate those eddy viscosity models available in the 
literature which were believed to offer the best chance of success in correlating the data. 
The eddy viscosity models chosen for this evaluation were those of Schetz (ref. 3) and 
Cohen (ref. 4) and that developed in a previously conducted air-air mixing study (ref. 6). 

These three eddy viscosity models were incorporated into the analysis of refer- 
ence 11. The validity of each model was determined by comparing analytical solutions to 
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hydrogen-air data from the present study and air-air mixing data from reference 6 . It 
was considered highly desirable that an eddy viscosity model permit the calculations to 
be initiated at the nozzle exit and to proceed continuously throughout the flow field. As 
eddy viscosity models are, in general, based upon dimensional analysis considerations, 
it was adjudged for this study that for the empirical constant to be truly a constant was 
too optimistic a viewpoint. Therefore, the value of the empirical constant associated with 
each model (assuming correct trends could be predicted) was chosen to best correlate the 
data. 

Data were obtained corresponding to the nonreactive turbulent mixing of a circular 
coaxial near -sonic jet of hydrogen surrounded by a supersonic parallel stream of air at 
M = 1.32 or M = 2.50. The hydrogen jet Mach number was approximately 0.9 for both 
airstream Mach numbers. Representative unit Reynolds numbers are 1.18 x 10 6 per 
meter for the hydrogen jet and 3.71 x 10 6 per meter and 1.39 x 10 7 per meter for the 
Mach 1.32 and 2.50 air jets, respectively. The jets mixed in an unconfined region at a 
static pressure of 1 atmosphere (1 atm = 1.013 x 105 N/m2). Both the hydrogen and air 
jets had total temperatures of approximately 300 K. Radial distributions of pitot pressure 
and hydrogen concentration were obtained downstream of the jet exit at various axial loca- 
tions. These data were reduced to velocity and hydrogen mass fraction profiles and are 
tabulated in appendixes A and B. 

SYMBOLS 

A area 

b a mixing-zone width defined as radial distance between points where veloc- 

ities are U 3 and U 4 

Cp local specific heat at constant pressure 

Cp a specific heat at constant pressure of pure air 

C n h ‘ specific heat at constant pressure of pure hydrogen 

d inner diameter of center nozzle, 11.6 mm 

f* parameter in Cohen's eddy viscosity model (see eq. ( 8 )) 

I ratio of integrated hydrogen flow rate to metered hydrogen flow rate (see 

eq. (2)) 
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empirical constant for use with the Z -difference eddy viscosity model (see 
eq. (3)) 


empirical constant for use with the unified eddy viscosity model (see eq. (6)) 

empirical constant for use with the Cohen eddy viscosity model (see eqs. (8) 
and (9)) 

empirical constant for use with the kinematic Z -difference eddy viscosity 
model (see eq. (13)) 

a characteristic length (see eq. (6)) 

local Mach number 

local molecular weight 

ratio of outer jet velocity to center-line velocity, Ua/uo 
metered center jet hydrogen flow rate 

integrated center jet hydrogen flow rate calculated from measured pressures, 
temperatures, and concentrations (see eq. (1)) 

velocity ratio fixed by turbulence level (see eq. (9)) 

turbulent Lewis number (the product of the turbulent diffusion coefficient and 
the constant-pressure specific heat divided by the turbulent thermal- 
conductivity coefficient) 

turbulent Prandtl number (the product of the constant-pressure specific heat 

and the turbulent viscosity divided by the turbulent thermal -conductivity 
coefficient) 

turbulent Schmidt number (ratio of turbulent Prandtl number to turbulent 
Lewis number) 

ratio of outer jet density to center -line density 



n l 


p 

R 

R 

T 


u 

U 1 

u 2 


u 3 


u 4 

x 


y 

z 


a 


P 


y 


€ t 


density ratio fixed by turbulence level (see eq. (9)) 

static pressure 

universal gas constant 

local gas constant 

static temperature 

total temperature 

axial velocity, m/sec 

velocity defined by equation (4) 

velocity defined by equation (5) 

velocity defined by equation (10) 

velocity defined by equation (11) 

axial coordinate 

radial coordinate 

a mixing-zone width defined as radial distance between points where velocities 
are u^ and U 2 

local mass fraction of hydrogen (local mass of hydrogen divided by sum of 
local mass of hydrogen and local mass of air) 

local volume fraction of hydrogen 

ratio of specific heats 

a displacement thickness defined by equation (7) 
eddy viscosity in kinematic form 
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dynamic eddy viscosity 


(Pe)t 

P density 

Subscripts: 

a evaluated in external airstream (see fig. 3) 

j evaluated in hydrogen flow at jet exit and on jet center line 

o evaluated on center line 

Bar over symbol indicates that parameter is nondimensionalized by the center 
nozzle external diameter (12.7 mm). 


APPARATUS 

Nozzle 

A sketch of the nozzle configuration used to generate the flow fields of this study is 
shown in figure 1(a). Nozzle contours are given in figures 1(b), 1(c), and 1(d). A Mach 
1.32 circular -contoured plug nozzle with a subsonic circular nozzle contained in the center 
body was used in the initial phase of this study. A second phase employed a Mach 2.50 
circular-contoured plug nozzle with the subsonic circular center nozzle. The exit diam- 
eter for both plug nozzles was 15.2 cm. The subsonic nozzles had inner and outer exit 
diameters of 11.6 mm and 12.7 mm, respectively. (See enlarged sketch in fig. 1(a)). The 
flow passage in the subsonic nozzles converged slightly in the flow direction to ensure that 
if choking occurred it would occur at the nozzle exit. The taper was approximately 
0.005 cm/cm on the diameter. The tapered section extended approximately 10 cm up- 
stream of the nozzle exit at which point the internal diameter was increased to 25.4 mm; 
this enlarged diameter resulted in the flow area being increased by a factor of approxi- 
mately 4. 

An exhaust duct 76 cm in diameter was located approximately 1 m downstream of 
the nozzle exit and thus provided means of exhausting the hydrogen-air jets from the test 
area. The inlet to the test cell was open to the atmosphere and provided sufficient flow 
area to prevent any detectable decrease in static pressure in the test chamber during data 
acquisition. 
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Survey Rake 


A remote controlled actuator was installed at various distances in the axial direc- 
tion from the nozzle exit and was used to traverse a survey rake across the flow field. 

The rake position was indicated by a servodriven counter. The counter which indicated 
probe location was calibrated to be 195 counts/cm with an observed uncertainty of 
±2 counts. The survey rake contained a pitot probe and a static-pressure probe which 
were 4.06 cm apart. Significant details of the rake and probes are shown in figure 2. A 
flow-field schematic is shown in figure 3. 

Gas Analysis 

Gas samples were piped directly from the rake to the injection valve in a gas chro- 
matograph. Standard gas chromatography techniques similar to those described in ref- 
erence 12 were used to perform the gas analysis. Briefly, this technique involves sepa- 
ration of the hydrogen from other gas components by passing the sample through a silica 
gel column followed by a molecular sieve column. The separated gas components then 
were passed through a thermistor detector. The output of this detector relative to -that 
of a reference detector exposed to pure nitrogen carrier gas is an indication of the volume 
of hydrogen in a given sample. The detector output was recorded on a strip chart 
recorder. Comparison of the peak height output to peak heights corresponding to known 
mixtures of hydrogen and nitrogen gave the volume percent hydrogen in a sample. Per- 
tinent data on the columns are as follows: 

Column 1: silica gel; mesh 60/70; diameter, 3.18 mm; length, 122 cm; conditioned for 
4 hr at 470 K with argon gas purge. 

Column 2: molecular sieve; 5 A; 70/80 mesh; diameter, 3.18 mm; length, 183 cm; 
conditioned for 6 hr at 620 K with argon purge. 

The nitrogen carrier gas flow rate was approximately 45 standard cm3/min. The columns 
were maintained at 310 K during gas analysis. Errors in gas analysis are not significant 
in comparison to the uncertainty of obtaining representative gas samples from the tur- 
bulent flow field. 


TEST PROCEDURE 

A constant flow of dry air at near-ambient temperature was supplied to the plug 
nozzle through a duct 36 cm in diameter. Ambient-temperature hydrogen was supplied 
to the inner nozzle through a support pipe. The temperature of each stream was mea- 
sured in the supply pipes by use of iron-constantan thermocouples. For the established 
flow conditions, static-pressure orifices near the outer nozzle exit indicated ratios of 
nozzle-exit static pressure to atmospheric static pressure of 1.02 and 0.97 for the Mach 
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1.32 and Mach 2.50 jets, respectively. Radial surveys of pitot pressure and hydrogen 
concentration were performed at several axial stations in the flow field. A limited num- 
ber of probe static-pressure surveys were also performed. Survey locations and tem- 
peratures for the two test conditions studied are summarized in table I. Schlieren photo- 
graphs were also taken to aid in interpreting the data. Typical flash and time-exposure 
schlieren photographs are presented in figures 4 and 5. 

During the early stages of the study, gas samples were extracted from the flow field 
through a pitot probe tip of the type shown in figure 2. Only the pitot pressure was used 
as the pumping force to extract these samples. In order to assess the accuracy of the 
data, the following integral was evaluated for each axial survey station: 


m-v 


= ^MdA 
\fwr 


(i) 


Details on the evaluation of equation (1) are given in appendix C. The magnitude of m x 
should be equal to the metered hydrogen flow rate m^, or the magnitude of I as given 
by 


I 



( 2 ) 


should be near unity if there are no substantial errors in the data. Any deviation of I 
from unity is due to experimental error. For hydrogen-air mixing, I is particularly 
sensitive to hydrogen concentration. Typical values of I computed from the initial data 
were in the range of 0.7 to 0,8 which indicated that 20 to 30 percent of the hydrogen jet 
flow was not accounted for in the gas sample data. Attempts to obtain higher concentra- 
tion values by aspirating the gas sample through the pitot probe and by pumping on the 
probe were unsuccessful inasmuch as no improvement was noted in the concentration data. 
A recent study (ref. 13) noted that differences between integrated and measured mass 
flows of 20 percent are considered typical. 

The problem of obtaining representative gas samples has been discussed in refer- 
ence 6. It was concluded in reference 6 that the actual physical mechanism which causes 
the sampling probe to obtain unrepresentative samples is not known; however, the results 
suggest that the erroneous concentration measurements are related to the local turbulence 
level in the flow field. It was noted in reference 14 that when sampling through a static 
probe, representative gas samples were always obtained. Therefore, several concentra- 
tion surveys were performed with the Mach 1.3 static probe (fig. 2) and a diaphragm -type 
pump to extract the gas sample. A comparison of the gas concentration obtained by pitot 
probe sampling and static-pressure probe sampling is presented in figure 6 for two axial 
stations. It is evident from the data in figure 6 that significantly higher gas concentrations 
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were obtained with the static probe sampling technique. Reduction of the static-pressure- 
probe sampling data of figure 6 resulted in ratios of integrated hydrogen flow rate to mea- 
sured hydrogen flow rate (eq. (2)) very near unity. Due to this excellent agreement, the 
static probe tips shown in figure 2 were used to extract all gas samples in the remainder 
of this test program. 

Some uncertainty in the axial location of the concentration data is introduced by 
sampling through the static probe. The uncertainty is due to the inability to positively 
determine whether the gas samples are representative of the flow at the probe tip or at 
the location of the static orifices. However, the boundary -layer flow on the probe will 
stabilize the small-scale turbulent flow field and thereby mixing is reduced in the vicinity 
of the probe surface. Thus, gas samples taken through the static orifices would be 
expected to have undergone less mixing than actually occurred in the undisturbed flow 
field. The concentration data are therefore expected to be representative of the flow 
field at some location ahead of the static orifice (but of course not ahead of the probe tip). 
For this study (including the data of fig. 6), the tip of the static probe was positioned at 
the same axial location as the tip of the pitot probe. The concentrations measured were 
assumed to be representative of the flow field and the axial location of the probe tip. A 
comparison of the dimensions of the static probes and the center nozzle diameter indi- 
cates that the concentration data may be displaced downstream from the pitot pressure 
survey by up to 1 jet diameter for the M a =1.32 data. Similarly, the concentration may 
be displaced up to 1.54 jet diameters from the pitot pressure survey for the Ma = 2.50 
data. The actual error introduced depends upon the axial total-pressure gradient and 
concentration gradient at each particular survey station and the effect of the interaction 
between the probe and the flow field. 

As an indication of the accuracy of the data, values of I from equation (2) have 
been listed in table I. As shown in table 1(a), values of integrated hydrogen flow rate 
ranged from 16 percent low to 4 percent high for the M a = 1.32 data. Poorer accuracy 
was achieved for the M a = 2.50 data where integrated hydrogen flow rates ranged from 
12 percent low to 29 percent high. It is noted that values of I greater than 100 percent 
were calculated in the region of the flow field where large gradients in pitot pressure and 
concentration existed (from x/d = 4.31 to x/d = 15.36). 

DATA REDUCTION 

The measured pitot pressures, total temperature, and volumetric concentration of 
hydrogen were reduced to velocity profiles and hydrogen mass fraction profiles. The 
method employed is identical to the technique used in evaluating equation (1). Specific 
equations are presented in appendix C. The center of the hydrogen concentration profiles 
was assumed to be the center of the flow field. An assumption of uniform static pressure 
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equal to atmospheric pressure was used in all data reduction. Probe static-pressure 
measurements were found to be seriously affected by shock-wave intersections with the 
probe and were not used in data reduction. 

It was observed during data reduction that the extent of the mixing zone (the region 
where hydrogen- air mixtures existed) could be determined from the fluctuations in the 
pitot-pressure signal. In the region of large gradients (in the mixing zone near the con- 
centration potential core), pitot-pressure fluctuations were as high as ±7.1 percent in an 
extreme case. (These fluctuations are not to be taken as indicative of the actual magni- 
tude of the turbulence of the flow as they are a function of the recording system employed, 
but the pitot-pressure fluctuations are an indication of large turbulent fluctuations in the 
mixing zone.) The fluctuation in pitot pressure decreased to an insignificant value as 
either the undisturbed hydrogen potential core flow or the undisturbed airflow was 
approached. Some evidence of large-scale mixing vortices, and thus associated fluctua- 
tions, may be seen in figures 4(a) and 5(a), as the mixing boundaries are irregular in the 
schlieren flash photographs. Mean values of pitot pressure from the strip chart record 
were used in data analysis. The local Mach number, mass fraction, and velocity data are 
tabulated in appendixes A and B. 

Concentration values required to compute local velocities were read from plots of 
the radial distribution of concentration inasmuch as pitot-pressure and concentration mea- 
surements were not necessarily made at the same radial location. 

THEORY 

Analysis 

The analysis of reference 11 was used to correlate the experimental data of this 
study and the air-air mixing data from reference 6. The analysis employs equilibrium 
chemistry, transformation techniques, and an explicit finite -difference technique to com- 
pute turbulent mixing and reacting of parallel streams of hydrogen and air. Axial pres- 
sure gradients and nonunity turbulent Prandtl and Lewis numbers are provided for in the 
analysis. The analysis is applicable to both the near and far field providing a proper 
eddy viscosity model is specified. Further details of the analysis are given in refer- 
ence 11. Information on the use of the computerized analysis may be found in reference 7. 

In order to analyze the hydrogen-air mixing data without allowing a reaction to 
occur, the airstream was simulated by pure nitrogen. For air-air mixing, the center jet 
was specified as an oxygen -nitrogen mixture corresponding to that of air (oxygen mass 
fraction 0.232, nitrogen mass fraction 0.768) and the outer stream as pure nitrogen. The 
oxygen thus defined the extent of mixing of the center jet. Specifying the concentration 
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composition as just noted retains the validity of the mixing analysis by maintaining the 
approximately correct molecular weight ratio between the streams. 

Eddy Viscosity Models 

Different eddy viscosity models were incorporated in the analyses to test their 
validity. The first model examined herein was developed in reference 6 and is expressed 
as 


(pe) t = kz(pu) Q (3) 

In equation (3), (pe) t is the dynamic eddy viscosity, k is an empirical constant, and 
(pu) Q is the mass flux per unit area on the jet center line. The mixing-zone width z 
is defined as the radial distance between the points where the local velocities are u^ 
and U£ as given by the following equations : 

u l = u a + 0.95^ - u a ) (4) 

u 2 = u a + 0.50^u o - u a ) (5) 

The definitions of uj, U 2 , and z are illustrated in the following sketch: 


u 
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The eddy viscosity model of equation (3) is hereinafter referred to as the Z -difference 
model and was satisfactorily used in reference 6 to correlate air-air mixing data. Val- 
ues of k employed in reference 6 were of the order of 0.01 and varied slightly with test 
conditions. 

The second model considered relates the eddy viscosity to the mass flow defect (or 
excess) across the mixing zone. This model was termed the "unified” eddy viscosity 
model in reference 3 and is given by 


(Pe) t = 


77 k 2 (pu)a^ 

L 


(6) 


In equation (6), ^kg is an empirical constant, (pu) a is the mass flow per unit area in 
the outer stream, and L is a characteristic length assumed to be the nozzle radius. 
The displacement thickness 6* is expressed as 



(7) 


Correlation of data using the unified model of equation (6) was reported in reference 3 by 
using a value of rrkg of 0.018. In the context of the usage of equation (6) in reference 3, 
the constant 7*2 was unchanged which implies that the constant is independent of test 
conditions and applicable to all data. A limitation in equation (6) is that defining the char- 
acteristic length L as the nozzle radius is a very poor approximation in the region near 
the nozzle exit. Therefore, the application of equation (6) is expected to be confined to 
the downstream region of the flow field. 

The third eddy viscosity model considered was developed in reference 4 and is 
hereinafter referred to as the Cohen viscosity model. It differs from the other two 
models by employing the kinematic eddy viscosity. The model is defined by the following 
equations: 


e * p o + p a ' 


, 0.8 


€ t =k3f ^r 


€ t = k 3 


1 + n l\°’ 8 . 


1 - mi 


\U r 


^1 + n^(l + mn) 
(1 + n) ^1 + rn-^n^ 


( 8 ) 

(9) 


In equations (8) and (9), kg is an empirical constant, f* is an empirical parameter 
equal to unity for incompressible mixing but which may vary with Mach number, p Q is 
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the density on the jet center line, p a is the density in the external flow, u Q is the 
velocity on the jet center line, and u a is the velocity in the external flow. In equa- 
tions (8) and (9), the mixing-zone width b is defined as the distance across the mixing 
zone between the points where the velocities are u 3 and U4 as given by the following 
equations: 


u 3 = u a + 0.95(u o - u a ) (10) 

u 4 = u a + 0.05^u o ~ u a) ^ ^ 

In equation (9) , n is the density ratio P a /P 0 > m is the velocity ratio u a ^/u 0 , nj is 
a density ratio fixed by turbulence level, and mi is a velocity ratio fixed by turbulence 
level. Equation (8) is to be used if m5mj,and equation (9) is to be used for m>m 1 . 
The value of is to be calculated at the axial station where m = mp If the initial 
velocity ratio exceeds mj, n^ is to be taken as the ratio of the external stream density 
to the initial jet density. In reference 4, f* was taken as unity and mj was taken equal 
to 0.40. A value of k 3 of 0.00764 was suggested for the core region and a somewhat 
larger value of 0.0089, for the downstream region. As it did not appear possible to change 
the value of the constant k 3 between the near field and the far field without introducing 
discontinuities, a constant value of k 3 was used throughout the flow field for all com- 
putations presented herein. Note that the dynamic eddy viscosity (pe) t , obtained from the 
local density and equations (8) and (9), varies in both the radial and axial directions. The 
other two eddy viscosity models as given by equations (3) and (6) permit eddy viscosity 
variations in only the axial direction. 


DATA PRESENTATION AND CORRELATION 
General Comments 

It is recognized that a transition region exists between the quasi -two -dimensional 
near-field mixing region and the far -field fully developed profiles. However, turbulent 
flow theory is not sufficiently developed to provide a means of treating this transition 
region. Therefore, in this study, as has been the approach of most investigators (see 
ref. 3 for further discussion on the neglect of the transition region), the flow field is con- 
sidered to consist of only a near-field mixing region and a far-field mixing region (fig. 3). 
The neglect and inadequate knowledge of how to treat the transition region result in poorer 
data correlation in the region near the end of potential core than in the remainder of the 
flow field, as evidenced by the results of reference 6 and of this study. This inaccuracy 
has been accepted in order to provide an eddy viscosity model which permits initiation of 
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calculations at the nozzle exit and which permits continuous calculations to proceed 
throughout the flow field. 

A primary purpose of this study was to evaluate several eddy viscosity models. It 
was realized that limiting the theoretical calculations to the two sets of hydrogen-air 
mixing data generated in this study was too restrictive an evaluation. Therefore, air -air 
mixing data from reference 6 was included in this study to evaluate the three eddy vis- 
cosity models of interest. 

A technique of evaluating a particular eddy viscosity model is to compare experi- 
mental center -line velocity distributions with those predicted by a particular model and a 
selected value of the empirical constant. Once the empirical constant is selected, it 
remains to determine the best values of the turbulent Lewis number Nj__ e and turbulent 
Prandtl number Np r to correlate the center -line mass fraction distribution. It is 
important to remember that any eddy viscosity model can be made to fit a particular 
center-line velocity data point by a judicious choice of the empirical constant. However, 
the remainder of the data distribution will not necessarily be correlated. It was found in 
a previous study (ref. 6) that if the center-line velocity and concentration axial distribu- 
tions were reasonably well correlated, then the radial profiles were also reasonably cor- 
related. Therefore, the center -line correlation technique was used to evaluate eddy vis- 
cosity models in this study. Center-line velocity data from the current hydrogen-air 
mixing study and from the air-air mixing study of reference 6 are first presented. Pres- 
entation of radial velocity profiles and center-line and radial distributions of hydrogen 
mass fraction from the hydrogen-air mixing data are deferred until the eddy viscosity 
models are evaluated. 


Hydrogen-Air Velocity Data 

M a = 1.32 .- The center -line velocity data for M a = 1.32 and Mj = 0.89 are pre- 
sented in figure 7(a). It is noted that the center-line velocity decays in a consistent man- 
ner and approaches the free-stream velocity at the most downstream station (x/d = 63.6). 
Integrated mass flows as previously discussed and given in table I give a good degree of 
confidence in the reliability of these data. It was noted during data reduction that the data 
for x/d = 0 were for a higher hydrogen jet total temperature than the remainder of the 
data. Therefore, the total temperature (and thus the velocity) in the hydrogen jet was 
reduced to a value representative of the entire data set before initiating theoretical 
calculations. 

M a = 2,50 .- The center-line velocity data for M a = 2.50 and Mj = 0.91 are pre- 
sented in figure 7(b). For x/d = 4.31, 8.75, and 15.36 integrated mass flows averaged 
27 percent high. (See table I.) It is believed that the high values resulted from uncertain 
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concentration measurements in the region of large radial gradients where large fluctua- 
tions occur. The accuracy of the data at the downstream stations (x/d = 19.8 to 58.0) is 
believed to be good inasmuch as integrated mass flows of 0.88 to 0.96 were calculated. 

The data show an unexpected decrease in center-line velocity to values below the 
f ree _stream velocity for values of x/d near 20 and larger. Calculations indicate that the 
concentration measured near the center line at x/d = 19.8 would have to be increased 
13.6 percent in order to increase the computed center-line velocity to the magnitude of 
the free-stream velocity. Uncertainties in the concentration measurements are not 
believed large enough to explain the below free-stream velocity values. Center-line 
velocity data at stations beyond x/d = 19.8 are less sensitive to the concentration mea- 
surements due to the low magnitude of the concentration. 

The unexpected low center -line velocities may be explained in terms of the relative 
rates of mass and momentum transfer. The fact that mass transfer is faster than 
momentum transfer has been established experimentally (turbulent Schmidt number less 
than 1 for similar profiles). Therefore, the velocity can decrease axially even though the 
Mach number is increasing, as was found for the Mach 1.32 hydrogen-air velocity distri- 
bution. If the mass transfer is rapid enough relative to the momentum transfer, then a 
decrease in center-line velocity below free-stream velocity can occur as exemplified by 
the Mach 2.50 data. The difference in trends between the Mach 2.50 and Mach 1.32 
center-line velocity data may be explained by noting that there was a smaller difference 
between jet and free-stream velocities for the Mach 2.50 data. Therefore, relatively less 
rapid momentum transfer would occur and the mass transfer effect would be more pre- 
dominant than in the Mach 1.32 hydrogen-air mixing data. 

A complex pattern of shock waves is evident in the schlieren photographs of figure 5. 
Inasmuch as no significant degradation in free-stream Mach number occurred over the 
entire survey length of 58 diameters, it is believed that the flow disturbances were small 
and did not significantly affect flow-field development. 

Air-Air Velocity Data 

M a = 1.268.- An air-air mixing study of circular, coaxial, parallel compressible air 

jets has been reported in reference 6. The jets had total temperatures of approximately 
300 K and exhausted to the atmosphere. Mixing of the jets was studied with the use of 
tracer gas. Figure 7(c) presents the center-line velocity data for a Mach 0.813 inner jet 
surrounded by a Mach 1.268 outer jet. Large initial boundary layers were prevalent in 
both jets and caused the initial decrease in center-line velocity shown in figure 7(c). 
Further details about these data and the air-air data in figure 7(d) may be found in ref- 
erence 6. 
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M a = 1.302.- Center-line velocity data from reference 6 for air-air mixing of a 

Mach 0.942 circular coaxial jet surrounded by a Mach 1.302 jet is given in figure 7(d). 
The test conditions differ from those for the air data in figure 7(c) only in the slightly 
larger values of M a and Mj and the use of a larger center nozzle. Data obtained at 
x/d = 49 are reported to be of uncertain accuracy because the outer mixing boundaries 
had extended into the mixing zone; this necessitated large corrections to the data at this 
station. 


Comparison of Unified and Z -Difference Viscosity Models 

Theoretical center-line velocity distributions for hydrogen-air mixing and air-air 
mixing just discussed were generated by using the unified and Z -difference viscosity 
models and are presented in figure 7. Computations were started at the nozzle exit in 
two cases in order to support the contention that the unified model would lose its validity 
near the exit of an axisymmetric nozzle. 

It is evident from inspection of figure 7(a) that the unified model drastically over- 
predicts the rate of decay of the center-line velocity for M a = 1.32 when computations 
are started at the nozzle exit plane. Computations were also started at a downstream 
station (x/d = 9.58) with the initial conditions obtained from the actual measured profiles. 
Again, the unified model predicts much too rapid mixing when the recommended constant 
7rk2 of 0.018 was used. Calculations indicated that a constant 7rk2 of 0.00306 was a 
more reasonable value for use with the unified model and the present data as shown by 
the dashed curve in figure 7(a). 

The measured profiles at x/d = 9.58 were also used to initiate computations using 
the Z -difference model for the hydrogen-air data of figure 7(a) and a constant k value 
of 0.028. It is evident that the Z -difference model is superior to the unified model for 
these data. 

Figure 7(b) presents the results of using the measured profiles at x/d = 4.31 to 
initiate calculations for the M a = 2.50 data with both the unified and Z -difference 
models. Again the unified model with its recommended constant ?7k2 of 0.018 greatly 
overpredicts the center-line velocity decay. Calculations indicated that a value of ?rk 2 
of 0.00101 was more applicable to the data as evident in the right-hand curve of fig- 
ure 7(b). 

The center-line velocity decay predicted by using the Z-difference model and a 
constant k of 0.0255 is also shown in figure 7(b). As for the Mach 1.32 data of fig- 
ure 7(a), the unified model prediction deviates more sharply from the data than the 
Z-difference model prediction. Neither model predicts the decay of center-line velocity 
below free-stream velocity as exhibited by the data. A discussion as to why the calcula- 
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tions do not predict the center-line velocities below free-stream velocity may be found 
in the section entitled "M a = 2.50 Hydrogen-Air Radial Profiles." 

The center-line velocity distribution predicted for the M a = 1.268 air data by 
using the unified model with calculations initiated at the nozzle exit is shown in fig- 
ure 7(c). As for the hydrogen-air data, the predicted center-line velocity change is much 
too rapid. Since the velocity trends appeared to be correct, the unified model solution 
was calculated for a constant 7rk 2 of 0.00356. Comparison of this solution with the 
Z-difference model solution (k = 0.0078) shown in figure 7(c) clearly shows the superiority 
of the Z-difference model. 

Theoretical center-line velocity distributions for the M a = 1.302 air data are pre- 
sented in figure 7(d). Calculations were initiated at x/d = 17.2 and used the experimen- 
tally determined profiles from reference 6. The unified model together with the recom- 
mended value of 7 rk 2 of 0.018 again overpredicts the mixing rate and associated velocity 
increase. (Similar calculations initiated at x/d = 25 are presented in fig. 12 of ref. 3. 

A similar trend of overpredicting the mixing rate was exhibited.) When the constant 7 rk 2 
was reduced to 0.0094, the unified model solution was only slightly less satisfactory than 
the Z-difference model solution with k = 0.0098, as evidenced by the curves in figure 7(d). 

For all data considered, it is concluded that the Z-difference model correlates data 
more satisfactorily than the unified model. The recommended value of the empirical con- 
stant 7 rk 2 of 0.018 for use with the unified model was unsatisfactory and too large by a 
factor of approximately 2 to 6. The unified model may find an application to the mixing 
of streams of nearly equal velocities across the mixing zone, where the mixing-zone width 
of the Z-difference model becomes indefinite, but where a mass flux difference exists 
between the streams. It is a desirable feature of any eddy viscosity model to be able to 
correlate data by initiating computations at the nozzle exit and to proceed downstream with 
the computations in a continuous manner. Satisfactory correlation and continuous com- 
putations from the nozzle exit were not possible with the unified model for the axisym- 
metric mixing data presented. Unsatisfactory correlation appears to result from an 
incorrect assumption that the viscosity is related to the mass flow defect or excess 
across the mixing zones. Continuous computations from the nozzle exit were not pos- 
sible because of a poor approximation to the characteristic length of the near-field mixing 
zone. 


Comparison of Z -Difference and Cohen Viscosity Models 

In figure 8, center -line velocity data of figure 7 are compared with theoretical 
center-line velocity distributions computed by using the Z-difference and Cohen viscosity 
models. All Cohen model solutions presented employ values of the constants f* and 
mj of 1.0 and 0.4, respectively, as recommended in reference 4. 
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The center-line velocity distribution corresponding to the M a = 1.32 hydrogen- 
air data and the Z -difference model is presented in figure 8(a). The computed velocity 
distribution is seen to be low in the region of jet exit, high in the immediate far field, and 
to only approach the data at the most downstream stations. 

The center-line velocity distribution computed by using the Cohen model and a value 
of IC3 of 0.002 is also presented in figure 8(a). (The values of turbulent Lewis number 
and turbulent Prandtl number quoted for the Cohen model calculations in figs. 8(a) and 8(b) 
are for later reference purposes.) Since the ratio of free-stream velocity to jet velocity 
was 0.366 (less than 0.4), equation (8) was used to initiate the Cohen model computations. 
As computations proceeded to an x/d of 6.81, the ratio of free-stream velocity to 
center-line velocity increased to 0.4; this necessitated a viscosity model change to equa- 
tion (9), as recommended in reference 4. The Cohen model is seen to give a much better 
correlation of the center-line velocity than the Z -difference model for the data of fig- 
ure 8(a). 

In figure 8(b), the result of computations using the Z -difference model is shown for 
the M a = 2.50 hydrogen-air data. The Z -difference model solution is seen to be low in 
correlating the velocity near the nozzle exit and high in the downstream region, as was 
the Z -difference solution shown in figure 8(a). 

The Cohen model solution, as computed by use of equation (9) for u a /uj > 0.4, is 
also shown in figure 8(b) for a value of k3 of 0.00129. The correlation with the 
M a = 2.50 hydrogen-air data is seen to be good. The solution deviates from the data 
only when the center-line velocity data decreases to values below the free-stream 
velocity. 

The Z -difference and Cohen model correlations for the air-air mixing data are 
shown in figures 8(c) and 8(d). The Cohen model solutions were computed by the use of 
equation (8) for u a /uj > 0.4. Both models correlate the air-air velocity data with approx- 
imately equal accuracy. Relative to each other, the Cohen model predicts slightly more 
rapid mixing in the near field and less rapid mixing in the far field than the Z -difference 
model. 

It is concluded from the discussion of the correlations in figure 8 that the kinematic 
viscosity model of Cohen is superior to the Z -difference model. However, the Cohen 
model required empirical constants ranging from 0.00129 to 0.0038 compared with a 
recommended value from reference 4 of approximately 0.008. Since the empirical con- 
stants were significantly different from the recommended value, a calculation was per- 
formed for conditions corresponding to a known solution of Cohen in reference 5. The 
results obtained from using the present calculation technique and that given in reference 5 
were found to be in good agreement, and the empirical constant was near 0.008. However, 
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the parameter mj (see eq. (9)) was taken to be 0.6 for this calculation, which has a 
direct effect on the value of the required constant. 

It was anticipated that the Z -difference model would correlate the hydrogen-air data 
in a more satisfactory manner than occurred, particularly since the model satisfactorily 
correlated the air-air data of reference 6. Inasmuch as the Cohen model was superior to 
the Z -difference model, a study was made to determine the significant differences in the 
two models. If one considers the Cohen model given by equation (9) and reduces it to 
primary variables, the result is the following expression: 
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In equation (12), the empirical constant has been combined with other parameters which 
are only a function of test conditions. As pointed out by Cohen in reference 5, the term 
in the brackets is not a strong function. Upon comparison with equation (3), it is noted 
that the two models differ primarily in the definition of the mixing width and in the use of 
a kinematic viscosity in the Cohen formulation. Calculations were made with the two dif- 
ferent mixing-width definitions. It was found that the mixing rate for a model was not 
significantly affected by changing mixing-width definitions, as long as the magnitude of 
the empirical constant was adjusted. It was concluded that the primary difference 
between the Cohen model of equation (9) and the Z -difference model was in the use of 
kinematic viscosity values in the Cohen model. Calculations were made to support this 
conclusion and the results are discussed in the next section. 


Comparison of Kinematic Z -Difference and Cohen Viscosity Models 
A kinematic form of the Z -difference model was developed and is given by 

et = k4zu 0 ( 13 ) 

In equation (13) all parameters are as previously defined with the exception of a new 
empirical constant k 4 . 

Figure 9 compares the results of computations using the kinematic Z -difference 
model and the Cohen viscosity model with the four sets of data previously discussed. In 
figure 9(a), the velocity distribution from the kinematic Z-difference model is seen to 
give excellent correlation of all data downstream of the velocity core for the M a = 1.32 
hydrogen-air data. A value of k 4 of 0.0175 was used in the correlation. The Cohen 
model correlation is seen to better predict the velocity core length but was less satisfac- 
tory than the kinematic Z-difference model in predicting downstream velocities. Com- 
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parison of the Z -difference model solution of figure 8(a) and the kinematic Z -difference 
model solution of figure 9(a) indicates that substantial better correlation was achieved 
with the kinematic Z -difference model. 

A comparison of the Cohen model solution of figure 8(a) for Np r = 0.8 and 
Nl 6 = 1.0 with the Cohen model solution of figure 9(a) for Np r = 1.0 and Nl g = 1.0 
shows only a very slight effect on the value of Np r . This comparison supports the pro- 
cedure of determining the empirical constant by correlating the velocity distribution and 
then determining the proper value of Np r and N^ e by correlating the concentration 
distribution. 

In figure 9(b), the velocity distribution resulting from computations using the kine- 
matic Z -difference model is seen to be slightly poorer than the Cohen viscosity model 
solution for the M a = 2.50 hydrogen-air data. Differences in trends between the Cohen 
and kinematic Z -difference model solutions are due almost entirely to the bracketed term 
in equation (12). Comparison of the Z -difference model solution of figure 8(b) and the 
kinematic Z -difference model solution of figure 9(b) again indicates the superiority of the 
kinematic Z -difference model. 

In figures 9(c) and 9(d), kinematic Z -difference model solutions and Cohen model 
solutions are presented for the air-air mixing data. For the air-air data, no significant 
difference between the kinematic Z -difference model solutions of figures 9(c) and 9(d) and 
the Z -difference model solutions of figures 8(c) and 8(d) is noted. The differences in the 
constants employed with the models (0.0078 as opposed to 0.0076 and 0.0098 as opposed to 
0.0106) are considered beyond the accuracy of the data and correlation method. 

On the basis of these comparisons, it is concluded that both the Cohen and the kine- 
matic Z -difference eddy viscosity models permitted correlation of the velocity distribu- 
tions satisfactorily and that a kinematic form of the eddy viscosity, which provides radial 
variations as well as axial variations in dynamic eddy viscosity through use of the local 
density, is superior to an eddy viscosity, which permits variations only in the axial direc- 
tion. The latter conclusion is based upon the fact that an eddy viscosity which employed 
the local density permitted correlation of both air-air and hydrogen-air mixing data, 
whereas an eddy viscosity which employed the center-line density satisfactorily corre- 
lated air-air mixing data only. The conclusion is in agreement with the results of ref- 
erence 15 which concluded that "the ultimate model for the turbulent transport coefficients 
must include variation in the axial as well as the radial direction.” The requirement for 
a kinematic eddy viscosity formulation found herein is also in agreement with a conclusion 
of reference 16 which reported a study of turbulent boundary layers with mass addition, 
combustion, and pressure gradients. In reference 16, it was concluded that agreement 
between experimental and predicted velocity profiles validated the assumption that the 
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Reynolds stress (which is directly related to the eddy viscosity by the velocity gradient) 
is kinematic in nature. 

Center-Line Mach Number and Hydrogen Mass Fraction Correlations 

Since both the kinematic Z -difference and the Cohen eddy viscosity models were 
satisfactory in correlating the velocity distributions, both models were used to correlate 
the center-line Mach number and the center-line hydrogen mass fraction distributions. 
These correlations, in particular the mass fraction distributions, determine representa- 
tive values of the turbulent Prandtl number and turbulent Lewis number. Concentration 
data corresponding to the air -air data need no further consideration inasmuch as the 
velocity correlations, and therefore the concentration distributions, are essentially iden- 
tical to those of reference 6. In reference 6, the concentration correlation for the air- 
air data was found to be insensitive to the actual magnitudes of Nj_, e and Np r and 
depended solely upon the turbulent Schmidt number Ng c . In reference 6 a turbulent 
Schmidt number of 0.6 was found to be representative of the M a = 1.268 and M a = 1.302 
data. 

The center-line Mach number distribution for M a = 1.32 hydrogen-air mixing is 
presented in figure 10(a). The center-line hydrogen mass fraction distribution for the 
same condition is presented in figure 10(b). Center-line Mach number and mass fraction 
distributions computed by using the Cohen and the kinematic Z-difference models are also 
shown in figure 10. In figure 10(a), the center -line Mach number distribution is best cor- 
related by the kinematic Z -difference model with NLe = 1-0 and Np r = 0.9 (Ng c = 0.9). 
The Cohen model solutions are seen to exhibit trends which diverge significantly from 
the data at x/d greater than approximately 15. Note that the two solutions shown for 
the Cohen model in figure 10(a) are for different values of Np r (0.8 and 1.0) but for val- 
ues of Nl 6 of unity. It is evident that the Cohen model solution with Nl 6 = 1.0 and 
N Pr = 1.0 best approximates the Mach data. 

In figure 10(b), the hydrogen mass fraction distribution is well correlated by the 
kinematic Z-difference model with Nj^ = 1.0 and Np*- = 0.9 for x/d greater than 
approximately 7. However, the Cohen model solution with Nl 6 = 1.0 and Np r = 1.0 
is seen to better predict the concentration core length. On the basis of the correlations 
of figures 9(a), 10(a), and 10(b), it is apparent that the kinematic Z-difference model cor- 
relates a larger portion of the flow field than the Cohen model. 

The center-line Mach number distribution for M a = 2.50 hydrogen-air mixing is 
presented in figure 10(c). The center-line hydrogen mass fraction distribution for the 
same test condition is presented in figure 10(d). It is noted in figure 10(c) that the kine- 
matic Z-difference model solution with N Le = 1.0 and Np r = 1.0 appears to best cor- 
relate the data. However, the kinematic Z-difference model solution with Nj^g = 1.0 
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and N Pr = 0.9, as was used in the correlations of figures 10(a) and 10(b), is believed to 
be within the accuracy of the data and correlation technique. In figure 10(c) the Cohen 
model solution with Nl g = 1.0 and Np r = 1.0 exhibits the same trend as in fig- 
ure 10(a), that is, diverging from the data in the region beyond x/d = 15. 

In figure 10(d), the mass fraction distribution is best correlated by the Cohen model 
solution with Nl 6 = 1.0 and Np r = 1.0. It is indeterminate whether values of 
N Le = N Pr = 1*0 o r N Le = 1.0, Np r = 0.9 are better for correlation using the 
kinematic Z -difference model mass fraction distributions of figure 10(d) because the 
N Le = 1*0, Np r = 1.0 solution results in better correlation of the mass fraction data 
for x/d less than approximately 13 and the N Le = 1.0, N Pr = 0.9 solution results in 
better correlation for data farther downstream. As in reference 6, the concentration cor- 
relation was found to be insensitive to the actual magnitude of N Le and Np r but 
depended only upon Ng c . The sole dependence upon Ng c is believed to be due to the 
fact that no significant transfer of heat is involved in the mixing process. 

In reference 8, values of N Sc ranging from 0.3 to 2.3 were reported for hydrogen- 
air, helium-air, and argon-air jets. (Some of the air jets considered in ref. 8 were 
heated.) No dependence of Ng c upon molecular weight was observed. Derivatives 
were obtained from experimental data in the determination of N gc in reference 8; thus, 
some scatter was undoubtedly introduced during computations due to data uncertainties. 
Reference 17 reported a value of Ng c near 0.7 for coaxial nearly equal density jets; 
this value closely agrees with results reported in reference 6. Reference 2 also reported 
a value of Ng c near 0.7 for carbon-dioxide — air and hydrogen-air jets in the downstream 
region. The various calculation techniques used by different investigators in determining 
N Sc make valid comparison difficult. As Ng c is a measure of the ratio of concentra- 
tion to velocity potential core length, it is difficult to conceive of N Sc as being other 
than near unity for hydrogen-air mixing. This follows from the fact that the velocity is 
extremely sensitive to hydrogen concentration and, thus, when the local hydrogen concen- 
tration changes, it would be closely accompanied by a change in velocity. 

M a = 1.32 Hydrogen-Air Radial Profiles 

Radial velocity profile data obtained from the M a = 1.32 hydrogen-air test are 
presented in figures 11(a) and 11(b). Radial hydrogen mass fraction profile data from 
the M a = 1.32 test are presented in figures 11(c) and 11(d). Correlations generated by 
using the Cohen viscosity model and the kinematic Z -difference model and the constants 
previously determined are also shown in the figures. 

The velocity data presented in figures 11(a) and 11(b) are seen to be generally self- 
consistent with the exception of the profile at x/d = 15.4 where some unexplained asym- 
metry is noted. The velocity is seen to be nearly uniform at x/d = 63.6; however, the 
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hydrogen mass fraction is approximately 1.6 percent on the center line indicating the total 
pressure distribution is not uniform. In figure 10(a) it is seen that the center-line Mach 
number at x/d = 63.6 is 1.18 as opposed to a free-stream Mach number of 1.32, which 
indicates mixing to a uniform condition has not been achieved. 

The velocity profile used to initiate the mixing calculations at x/d = 0.0 is shown 
as a solid line in figure 11(a). Because the data for x/d = 0.0 were for a higher tem- 
perature (and thus a higher velocity) than the remainder of the data, a lower more repre- 
sentative velocity profile was used to initiate the calculations. Analytical solutions were 
not possible beyond approximately 32 jet diameters from the nozzle exit because the 
mixing width used in the viscosity models fails as the velocity becomes uniform. There- 
fore, no correlations are presented for x/d = 42.8 and 63.6. With the exception of 
x/d = 5.51, the kinematic Z-difference viscosity model (dashed lines in figs. 11(a) and 
11(b)) is seen to give slightly better correlation than the Cohen viscosity model (solid 
lines in figs. 11(a) and 11(b)) for the M a = 1.32 velocity data. 

The hydrogen mass fraction profile data in figures 11(c) and 11(d) exhibit uniform 
trends and good repeatability throughout the range of the measurements. (Note the scale 
change between figs. 11(c) and 11(d).) The poorest data appear to exist near the profile 
maximum at x/d = 9.58 which is in the region of large axial gradients and large fluctua- 
tions immediately downstream of the potential core. 

As was found with the velocity correlation, the kinematic Z-difference solution 
results in slightly better correlation than the Cohen solution for all stations other than 
x/d = 5.51. 


M a = 2.50 Hydrogen- Air Radial Profiles 

Radial velocity profile data obtained from the M a = 2.50 hydrogen-air test are 
presented in figures 12(a) and 12(b). Radial hydrogen mass fraction profile data are pre- 
sented in figures 12(c) and 12(d). Correlations generated by using the kinematic 
Z-difference and Cohen viscosity models and the constants previously determined are 
also shown in the figures. 

The velocity profile for x/d = 0.0 (fig. 12(a)) is noted to have a slightly larger 
boundary layer in the airstream than the M a = 1.32 data (fig. 11(a)). In order to limit 
the number of input points to the computer program, and thus limit the computational 
time to a reasonable amount, the actual measured boundary layer was approximated by 
the solid line shown in figure 12(a). It was found during computations that a step velocity 
profile could be satisfactorily used as input as long as only center-line correlation was 
attempted. However, when considering correlation of radial profiles, the actual measured 
profile must be used as input. Some possible effect of the approximation of the initial 
profile may be seen in figure 12(a) for x/d = 4.31 and 8.75. However, the effect of the 
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wake resulting from the finite nozzle lip thickness is impossible to separate completely 
from the initial boundary -layer approximation. 

As was found during the center-line velocity correlation, the Cohen viscosity model 
solution is a better approximation of the M a = 2.50 hydrogen-air data. However, neither 
model predicts the reduction in center-line velocity to below free-stream velocity shown 
in figure 12(b) for x/d = 19.8, 37.3, and 58.0. A study of the finite-difference technique 
employed in the analysis (ref. 11) indicates that once a uniform velocity profile is 
obtained, no further change in velocity is permitted. Furthermore, the viscosity models 
fail as the velocity profile approaches uniformity, due to employing a mixing width based 
on velocity. Thus, both the analysis employed and the eddy viscosity models prevent any 
computation of center-line velocity values below free-stream values. 

The hydrogen mass fraction data in figures 12(c) and 12(d) exhibit wider spreading 
than the solutions predicted by either viscosity model, particularly for x/d = 8.75 and 
15.4. It is interesting that these two axial stations were high in integrated-mass-flow 
ratio. The theoretical solutions of figures 12(c) and 12(d) indicate a possible bias in the 
data toward high values of concentration in this region of large radial gradients. 

The Cohen viscosity model solution is seen to give better correlation than the kine- 
matic Z -difference viscosity model solution for the M a = 2.50 hydrogen-air data. Solu- 
tions were not possible for either viscosity model for x/d = 37.3 and 58.0 because the 
mixing width used in the models failed near x/d = 30 as the velocity profiles approached 
the uniform velocity condition. 

The data and empirical constants employed herein are summarized in table H. 
Solutions beyond the region of application of the two viscosity models used herein can be 
generated by use of a constant viscosity, inasmuch as the regions are far downstream of 
the nozzle exit and should be approaching the boundary condition of uniform viscosity at 
infinity. Such computations were not made herein because they are in a region where 
application to scramjet combustor design would involve interaction between adjacent jets 
and, thus, involve an entirely different physical problem. 

CONCLUDING REMARKS 

An investigation of the compressible turbulent mixing of coaxial concentric 
hydrogen-air jets has been conducted. Hydrogen mass fraction profiles and velocity pro- 
files were acquired for airstream Mach numbers of 1.32 and 2.50 corresponding to hydro- 
gen jet Mach numbers of 0.89 and 0.91, respectively. These two sets of hydrogen-air data 
and two sets of air-air data (the air-air data from a previous study) have been satisfacto- 
rily correlated. For the data correlated, the ratio of jet velocity to free-stream velocity 
ranged from 0.648 to 2.730 and the ratio of jet mass flux to free-stream mass flux per 
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unit area ranged from 0.647 to 0.0725. The total -temperature ratios of the streams for 
all data considered was near unity. 

Three different eddy viscosity models were incorporated into a finite-difference- 
type analysis to test their validity in correlating the data. The kinematic eddy viscosity 
model of Cohen and a modified form of an eddy viscosity model used in a previous study 
both satisfactorily correlated the data. 

A formulation of eddy viscosity which was previously employed to correlate air-air 
data and which only allowed axial variation in eddy viscosity was unsatisfactory in cor- 
relating the hydrogen-air data. It was found that a kinematic form of eddy viscosity which 
allowed radial variation as well as axial variation in dynamic eddy viscosity through 
incorporation of the local density was essential in order to correlate both the air -air and 
hydrogen-air data. 

A formulation of eddy viscosity developed by Schetz, termed the "unified eddy vis- 
cosity model," which is based on a mass flow defect (or excess) across the mixing region, 
proved unsatisfactory in correlating the axisymmetric data considered herein. 

Turbulent Schmidt numbers N Sc used to correlate the hydrogen-air data were of 
the order of 0.9 to 1.0. The air-air data have been previously correlated with a value of 
Ng c of 0.6. 

Empirical constants used to correlate the data are, as yet, an unknown function of 
test conditions. The empirical constants used with the eddy viscosity models used herein 
varied by a factor of between 2 and 3 for the data considered. Correlation of a large 
quantity of data should permit a determination of the proper value of the constant to 
employ for given test conditions. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., August 24, 1971. 
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APPENDIX A 


VELOCITY AND HYDROGEN MASS FRACTION FOR M a = 1.32 and M; = 0.89 a 

d J 


x/d = 5,51 x/d = 9.58 


y/d 

M 

u 

y/d 

M 

u 


y/d 

M 

u 

y/d 

a 

- 

y/d 

M 

u 

y/d 

a 

-6.595 

0.368 

126 

0.525 

0.374 

121 


-5.973 

1.330 

397 

-1.200 

o.ooc 


-3.587 

1.304 

394 

-1.584 

0.000 

-6.565 

.897 

289 

.534 

.823 

261 


-5.537 

1.330 

397 

-1.129 

.00( 


-3.137 

1.302 

394 

-1.381 

.000 

-6.538 

1.055 

331 

.543 

.996 

31{ 


-5.316 

1.328 

396 

-1.129 

.00( 


-2.709 

1.301 

393 

-1.235 

.001 

-6.476 

1.218 

371 

.565 

1.046 

32 f 


-5.126 

1.327 

396 

-.887 

.004 


-2.316 

1.304 

394 

-1.068 

.006 

-6.437 

1.284 

386 

.604 

1.138 

353 


-4.897 

1.327 

396 

-.816 

.014 


-1.813 

1.310 

396 

-1.068 

.006 

-6.317 

1.318 

394 

.631 

1.182 

363 


-4.429 

1.328 

396 

-.710 

.04C 


-1.394 

1.274 

387 

-.926 

.020 

-6.203 

1.324 

395 

.675 

1.227 

373 


-3.970 

1.337 

398 

-.635 

.079 


-1.332 

1.339 

402 

-.926 

.019 

-5.965 

1.326 

395 

.724 

1.251 

37£ 


-3.485 

1.338 

398 

-.547 

.146 


-1.187 

1.321 

402 

-.807 

.046 

-5.612 

1.324 

395 

.803 

1.280 

385 


-3.141 

1.338 

398 

-.547 

.143 


-1.125 

1.297 

401 

-.807 

.045 

-5.193 

1.322 

394 

.949 

1.297 

389 


-2.682 

1.339 

399 

-.454 

.223 


-1.107 

1.295 

402 

-.706 

.081 

-4.729 

1.322 

394 

1.050 

1.303 

390 


-2.237 

1.335 

398 

-.313 

.455 


-1.024 

1.247 

401 

-.706 

.082 

-4.262 

1.322 

394 

1.174 

1.308 

391 


-1.738 

1.334 

398 

-.141 

.792 


-1.001 

1.232 

401 

-.596 

.123 

-3.790 

1.317 

394 

1.381 

1.315 

393 


-1.354 

1.334 

398 

-.Ml 

.773 


-.922 

1.165 

404 

-.596 

.120 

-3.357 

1.328 

396 

1.654 

1.320 

394 


-1.103 

1.324 

396 

.000 

1.000 


-.869 

1.103 

408 

-.476 

. 177 

-2.881 

1.326 

395 

2.025 

1.322 

394 


-1.010 

1.297 

392 

,000 

1.000 


-.825 

1.044 

411 

-.476 

.182 

-2.471 

1.324 

395 

2.462 

1.326 

395 


-.904 

1.262 

389 

.119 

.873 


-.794 

.991 

408 

-.335 

.263 

-2.029 

1.320 

394 

2.956 

1.328 

396 


-.851 

1.220 

394 

.119 

.873 


-.701 

.902 

421 

-.335 

.261 

-1.628 

1.320 

394 

3.428 

1.328 

396 


-.825 

1.187 

404 

.194 

.731 


-.596 

.831 

436 

-.176 

.480 

-1.244 

1.309 

392 

3.829 

1.326 

395 


-.768 

1.104 

396 

.194 

.692 


-.503 

.782 

450 

-.035 

.504 

-.975 

1.294 

388 

4.301 

1.324 

395 


-.724 

1.035 

398 

.251 

.585 


-.335 

.758 

535 

-.035 

.480 

-.843 

1.271 

383 

4.835 

1.324 

395 


-.684 

.948 

395 

.432 

.218 


-.203 

.777 

677 

.097 

.492 

-.776 

1.245 

377 

5.157 

1.324 

395 


-.635 

.864 

400 

.432 

.224 


-.097 

.809 

726 

.097 

.492 

-.697 

1.204 

368 

5.629 

1.326 

395 


-.560 

.783 

425 

.525 

.109 


.057 

.822 

739 

.229 

.413 

-.649 

1.153 

355 

6.079 

1.327 

396 


-.441 

.755 

514 

.525 

.109 


.199 

.796 

674 

.229 

.413 

-.591 

1.079 

337 

6.481 

1.305 

391 


-.366 

.783 

616 

.578 

.063 


.322 

.782 

590 

.393 

.256 

-.565 

1.023 

323 

6.525 

1.196 

366 


-.313 

.810 

709 

.635 

.033 


.463 

.822 

509 

.534 

.140 

-.543 

.961 

306 

6.587 

.968 

308 


-.247 

.840 

807 

.666 

.019 


.534 

.891 

488 

.534 

.144 

-.538 

.823 

268 

6.617 

.520 

176 


-.384 

.779 

594 

.666 

.019 


.613 

.997 

475 

.613 

.091 

-.490 

.621 

781 





-.331 

.802 

680 

.715 

.008 


.644 

1.044 

463 

.693 

.049 

-.472 

.735 

911 





-.163 

.860 

927 

.816 

.001 


.701 

1.131 

449 

.693 

.049 

-.459 

.772 

952 





-.110 

.864 

1012 

.816 

.001 


.741 

1.209 

446 

.746 

.026 

-.432 

.817 

1001 





.009 

.860 

1060 

.882 

.000 


.807 

1.272 

430 

.754 

.030 

-.384 

.852 

1039 





.044 

.856 

1056 

.882 

.000 


.904 

1.323 

412 

.754 | 

.031 

-.340 

.862 

1049 





.124 

.840 

975 

.966 

.000 


1.050 

1.343 

404 

.843 

.012 

-.274 

.872 

1060 





.194 

.802 

858 




1.275 

1.345 

403 

.843 [ 

.014 

-.199 

.878 

1066 





.221 

.793 

807 




1.544 

1.345 

403 

.926 

.004 

-.154 

.882 

1070 





.265 

.779 

743 




1.937 

1.343 

403 

.926 

.004 

-.097 

.884 

1072 





.313 

.774 

676 




2.418 

1.339 

402 

1.063 

.000 

-.018 

.886 

1074 





.344 

.774 

632 







1.063 

.000 

.119 

.886 

1074 





.424 

.815 

543 







1.187 

.000 

.251 

.876 

1064 





.481 

.927 

528 









.322 

.868 

1056 





.534 

1.087 

520 









.393 

.846 

1032 





.587 

1.184 

485 









.463 

.760 

938 





.618 

1.227 

466 









.485 

.615 

773 





.688 

1.303 

431 
















.812 

1.332 

400 
















1.090 

1.338 

398 
















1.518 

1.340 

399 
















1.998 

1.343 

400 
















2.515 

1.345 

400 
















2.973 

1.344 

400 












i 


— - 



. 

3.428 

1.343 

400 










a Total temperatures of the air and hydrogen jets for each x/d are given in table 1(a). 

^Velocities tabulated for x/d = 0.0 were computed for a hydrogen jet total temperature of 295 K, inasmuch as this 
value was more representative of the entire data rather than the actual measured value of 306 K. 
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APPENDIX A - Concluded 


x/d = 15.44 


y/d 

M 

u 


y/d 

a 

4.253 

1.331 

398 


- 1.504 

0.000 

4.068 

1.330 

397 


- 1.068 

.012 

3.768 

1.325 

396 


- 1.068 

.013 

3.335 

1.325 

396 


- 1.068 

.012 

2.991 

1.331 

398 


-.878 

.031 

2.550 

1.331 

398 


-.878 

.031 

2.268 

1.330 

397 


-.613 

.084 

1.871 

1.327 

397 


-.613 

.084 

1.513 

1.319 

396 


-.494 

.114 

1.328 

1.302 

397 


-.362 

.154 

1.213 

1.269 

398 


-.216 

.204 

1.103 

1.203 

394 


-.009 

.235 

-.997 

1.130 

393 


.150 

.225 

-.887 

1.053 

394 


.309 

.181 

-.807 

.988 

393 


.543 

.106 

-.688 

.918 

409 


.741 

.056 

-.543 

.851 

427 


.988 

.015 

-.441 

.826 

447 


.988 

.015 

-.322 

.820 

479 


1.116 

.005 

-.150 

.820 

528 


1.116 

.005 

.075 

.830 

551 


1.209 

.002 

.256 

.842 

529 


1.209 

.002 

.410 ' 
.525 
.631 
.715 
.803 
.887 
1.019 
1. 271 
1.509 
1.932 
2.400 
2.846 

.893 

.959 

1.066 

1.141 

1.218 

1.269 

1.311 

1.327 

1.330 

1.331 
1.331 
1.331 

499 
483 
483 
478 
1 467 
«> 
424 
! 399 

397 

398 
398 
398 


1.407 

.000 


x/d = 25.2 


y/d 

M 

u 


y/d 

a 

- 4.570 

1.313 

389 


- 2.444 

0.000 

- 4.068 

1.322 

391 


- 2.021 

.000 

- 3.573 

1.325 

392 


- 1.593 

# .001 

- 3.176 

1.326 

392 


- 1.390 

.005 

- 2.797 

1.327 

392 


- 1.390 

.005 

- 2.303 

1.327 

392 


- 1. 156 

.011 

- 1.862 

1.323 

392 


-.940 

.026 

- 1.372 

1.294 

397 


-.940 

.026 

- 1.200 

1.259 

401 


-.715 

.045 

- 1.072 

1.216 

404 


-.715 

.045 

-.957 

1.172 

411 


-.512 

.066 

-.891 

1.142 

412 


-.512 

.066 

-.785 

1.082 

410 


-.234 

.090 

-.657 

1.032 

415 


-.234 

.089 

-.547 

.996 

421 


-.071 

.102 

-.371 

.947 

428 


-.071 

.100 

-.097 

.919 

443 


.141 

.098 

.137 

.922 

443 


.141 

.099 

.270 

.933 

438 


.371 

. 0B5 

.446 

.960 

431 


.371 

.086 

.600 

1.005 

422 


.587 

.062 

.719 

1.049 

414 


.587 

.062 

.821 

1.126 

421 


.812 

.038 

.957 

1.193 

417 


.812 

.038 

1.081 

1.241 

411 


1.072 

.017 

1.209 

1.281 1 

404 

1 

1.072 

.016 

1.407 

1.311 

397 

| 

1.244 

.007 

1.619 

1.324 

395 

i 

1.526 

.002 

2.007 

1.327 1 

393 


1.526 

.002 

2.528 

1.328 1 

393 

1 

1.932 

.000 

2.943 

1.327 ! 

392 

1 

1.932 

.000 

3.335 

1.325 | 

392 




3.953 

1.321 

391 

1 




x/d = 42.8 


y/d 

M 

u 


y/d 

a 

- 7.350 

0.728 

240 


- 3,282 

0.000 

- 7.266 

.829 

270 


- 3.282 

.000 

- 6.432 

.949 

303 

• 

- 2.876 

.000 

- 6.088 

1.038 

327 


- 2.876 

.000 

- 5.625 

1.143 

353 


- 2.413 

.000 

- 5.409 

1.176 

362 


- 2.413 

.000 

- 5.170 

1.223 

373 


- 2.021 

.001 

- 4.853 

1.265 

383 

,- 2.021 

.001 

- 4.363 

1.306 

392 

; - 1.584 

.006 

- 3.860 

1.324 

396 

1 - 1.584 

.005 

- 3.485 

1.332 

398 


- 1.116 

.016 

- 3.181 

1,334 

398 


- 1.116 

.015 

- 2.793 

1.334 

398 


-.671 

.030 

- 2.290 

1.330 

398 


-.671 

.029 

- 1.884 

1.314 

399 


-.190 

.041 

- 1.544 

1.278 

401 


-.190 

.041 

- 1.249 

1.223 

401 


-.031 

.041 

- 1.010 

1.164 

399 


-.031 

.041 

-.754 

1.114 

401 


.216 

.040 

-.326 

1.073 

412 


.657 

.028 

-.040 

1.068 

416 


1.125 

.016 

.212 

1.087 

418 


1.535 

.006 

.428 

1108 | 

415 


1.557 

.005 

.710 

1.181 

421 


1.968 

.001 

.838 

1.217 

423 


1.968 

.001 

1.063 

1.257 

421 


2.413 

.000 

1.288 

1.293 

416 


2.413 

.000 

1.654 

1.324 

406 


2.828 

.000 

2.082 

1.332 

400 




2.515 

1.332 

398 




3.031 

1.330 

397 




3.432 

1.320 

395 




3.957 

1.299 

390 




4.443 

1.259 

381 




4.879 

1.210 

370 

1 




x/d = 63.8 


y/d 

M 

u 


y/d 

a 

- 8.091 

0.633 

212 


- 4.244 

0.000 

- 7.650 

,698 

232 


- 3.763 

.000 

- 7,187 

.769 

253 


- 3.763 

.000 

- 6.776 

.827 

270 


- 3.388 

.000 

- 6.304 

.898 

290 


- 2.925 

.001 

- 5.872 

.967 

309 


- 2.501 

.002 

- 5.457 

1.037 

327 


- 2.043 

.004 

- 5.065 

1.096 

342 


- 1.566 

.007 

- 4.690 

1.145 

354 


- 1.116 

.010 

- 4.293 

1.193 

366 


-.675 

.014 

- 3.851 

1.250 

380 


-.675 

.015 

- 3.582 

1.270 

384 


-.238 

.017 

- 3.221 

1.285 

388 


-.238 

.017 

- 2.929 

1.291 

391 


.172 

.017 

- 2.607 

1.296 

393 


.618 

.015 

- 2.188 

1.291 

396 


1.046 

.011 

- 1.760 

1.281 

400 


1.500 

.008 

- 1.434 

1.259 

400 


1.954 

.004 

- 1.129 

1.237 

401 


2.404 

.002 

-.874 

1.214 

400 


2.837 

.001 

-.631 

1.199 

401 


3.265 

.000 

-.371 

1.188 

402 


3.710 

.000 

-.106 

1.181 

402 


4.165 

.000 

.154 

1.186 

402 


4.597 

.000 

.379 

1. 186 | 

400 




.622 

1.199 

401 




.878 

1.219 

402 




1.121 

1 . 24 ! 

403 




1.354 

1.257 

402 




1.654 

1.274 

400 




1.915 

1.285 

398 




2.294 

1.291 

395 




2.612 

1.287 

391 




3.018 

1.274 

387 




3.463 

1.255 

381 




3.900 

1.224 

374 




4.363 

1.179 

363 




4.672 

1.133 

351 ! 




5. 166 

1.070 

336 | 




5.598 

.998 

317 




6.075 

.923 

297 ; 




6.604 

.843 

274 




7.147 

.764 

251 




7.747 

.658 

219 




8.365 

.587 

197 j 






APPENDIX B 


VELOCITY AND HYDROGEN MASS FRACTION FOR M a = 2.50 AND Mj = 0.91 a 

, X/d ° 0,0 x/dr4 31 x/d - 8.75 x/d - 15.36 


y/d 

M 

u 

y/d 

M 

u 


y/d 

M 

u 

y/d 

a 

- 6.666 

0.432 

153 

0 . 45 < 

0 . 59 ? 

76 


- 2.868 

2.588 

604 

- 0 . 97E 

0.000 

- 6.591 

.932 

311 

.472 

.432 

551 


- 2.673 

2.567 

602 

-. 86C 

.001 

- 6.512 

1.548 

461 

.485 

. 30E 

40 : 


- 2.550 

2.528 

598 

-.768 

.005 

- 6.437 

2.105 

555 

.so: 

1.096 

35 F 


- 2.334 

2.503 

596 

-.706 

.019 

- 6.159 

2.499 

603 

.516 

1.278 

40 


- 2.201 

2.536 

599 

-.662 

.046 

- 6.026 

2.502 

603 

.543 

1.522 

455 


- 2.012 

2.593 

605 

-.587 

.114 

- 5.572 

2.494 

603 

.569 

1.673 

483 


- 1.831 

2.598 

605 

-.587 

.115 

- 5.325 

2.531 

606 

.596 

1.760 

50C 


- 1.610 

2.577 

603 

-.529 

.197 

- 4.862 

2.538 

607 

.622 

1.834 

513 


- 1.518 

2.562 

602 

-.529 

.194 

- 4.412 

2.531 

606 

.688 

1.975 

536 


- 1.407 

2.531 

599 

-.472 

.316 

- 3.922 

2.538 

607 

.768 

2.092 

553 


- 1.310 

2.489 

594 

-.472 

.316 

- 3.468 

2.565 

610 

.860 

2.204 

568 


- 1.218 

2.414 

586 

-.375 

.549 

- 3.168 

2.586 

612 

.962 

2.287 

579 


- 1.147 

2.307 

574 

-.300 

.773 

- 2.943 

2.511 

604 

1.054 

2.333 

584 


- 1.072 

2.214 

562 

-.194 

.959 

- 2.493 

2.520 

605 

1.147 

2.379 

590 


-.949 

2.135 

552 

-.079 

1.000 

- 2.215 

2.545 

608 

1.253 

2.434 

596 


-.874 

2.066 

543 

.066 

1.000 

- 1.990 

2.524 

606 

1.368 

2.472 

600 


-.785 

1.967 

540 

.066 

1.000 

- 1.562 

2.503 

603 

1.522 

2.507 

604 


-.701 

1.835 

569 

.150 

.945 

- 1.394 

2.461 

599 

1.703 

2.527 

606 


-.649 

1.675 

625 

.234 

.841 

- 1.284 

2.423 

595 

1.932 

2.531 

606 


-.596 

1.424 

656 

.313 

.692 

- 1.134 

2.363 

588 

2.126 

2.519 

605 


-.560 

1.206 

642 

.424 

.402 

-.904 

2.240 

573 

2.303 

2.496 

603 


-.521 

1.025 

639 

.503 

.255 

-.746 

2.067 

549 

2.625 

2.481 

601 


-.476 

.922 

671 

.551 

.168 

-.578 

1.732 

496 

2.912 

2.472 

600 


-.397 

.874 

778 

.591 

.115 

-.547 

1.600 

471 

3.238 

2.575 

611 


-.326 

.889 

921 

.649 

.053 

-.499 

.718 

248 

3.684 

2.541 

607 


-.229 

.897 

1051 

.693 

.026 

-.405 

.309 

403 

4.147 

2.506 

604 


-.053 

.897 

1092 

.715 

.017 

-.468 

.507 

651 

4.663 

2.514 

605 


■ .106 

.897 

1092 

.759 

.007 

-.454 

.644 

815 

5.237 

2.493 

602 


.212 

.897 

1092 

.860 

.001 

-.446 

.662 

836 





.375 

.866 

797 

.860 

.001 

-.432 

.718 

901 





.468 

.911 

676 

.944 

.000 

-.415 

.783 

973 





.521 

1.040 

656 



-.397 

.819 

1013 





.534 

1.081 

651 



-.379 

.857 

1054 





.551 

1.182 

658 



-.344 

.877 

; 1075 





.574 

1.280 

657 



-.300 

.896 

1095 





.626 

1.534 

628 



-.053 

.907 

1107 





.657 

1.715 

609 



.154 

.907 

1107 





.693 

1.805 

577 



.251 

.907 

1107 





.719 

1.871 

560 



.353 

.871 

1068 





.776 

1.947 

542 



.366 

.865 

1062 





.851 

2.040 

542 



| .384 

.845 

1040 





.931 

2.097 

546 



J .406 

.815 

1008 





1.001 

2.160 

555 



.424 

.769 

957 





1.076 

2.202 

560 



.446 

.668 

843 





1.178 

2.303 

573 










1.266 

2.446 

590 










1.407 

2.542 

600 










1.548 

2.579 

604 










1.840 

2.588 

604 










2.109 

2.562 

602 










2.409 

2.484 

594 










3.044 

2.547 

600 










3.335 

2.525 

598 










3.706 : 

2.557 

601 









_ 

4.010 : 

2.557 

60lJ 




a Total temperatures of the air and hydrogen jets for each x/d are given In table f(b). 


y/d 

M 

u 

y/d 

a 


y/d 

M 

u 

y/d 

Of 

-4.173 

2.52E 

590 

- 1.07( 

O.OOC 


-4.27 

2.515 

592 

-1.429 

0.000 

-3.865 

2.560 

594 

-.776 

.015 


-3.60 

2.53 

593 

-1.204 

.004 

-3.591 

2.555 

593 

-.67E 

.058 


-2.98 

2.515 

592 

-1.204 

.004 

-3.018 

2.545 

592 

-.543 

.160 


-2.60 

2.51C 

591 

' -.984 

.009 

-2.722 

2.545 

592 

-.543 

.165 


-2.315 

2.503 

590 

-.768 

.028 

-2.594 

2.616 

599 

-.419 

.298 


-1.90 

2.489 

589 

-.657 

.063 

-2.400 

2.591 

597 

-.419 

.300 


-1.575 

2.441 

584 

-.657 

.064 

-2.113 

2.555 

593 

-.247 

.500 


-1.376 

2.386 

583 

-.525 

.103 

-1.641 

2.518 

589 

-.101 

.651 


-1.25" 

2.336 

585 

-.331 

.198 

-1.456 

2.487 

586 

.079 

.591 


-1.101 

2.255 

586 

-.106 

.252 

-1.328 

2.407 

577 

.296 

.44! 


-1.050 

2.190 

582 

.110 

.240 

-1.156 

2.329 

569 

.476 

.256 


-.957 

2.097 

582 

.331 

.196 

-1.041 

2.237 

560 

.728 

.035 


-.847 

1.912 

581 

.556 

.106 

-.909 

2.148 

561 

.869 

.006 


-.803 

1.783 

574 

.657 

.058 

-.856 

2.066 

558 

.979 

.001 


-.706 

1.595 

597 

.759 

.034 

-.807 

1.988 

559 

1.103 

.000 


-.649 

1.377 

576 

.882 

.011 

-.750 

1.800 

573 




-.521 

1.206 

603 

1.010 

.004 

-.706 

1.639 

590 




-.406 

1.067 

612 

1.107 

.001 

-.688 

1.503 

576 




-.159 

.995 

658 

1.209 

.000 

-.600 

1.179 

584 




.013 

1.001 

656 



-.534 

1.023 

595 




.176 

1.051 

658 



-.300 

.918 

776 




.388 

1.170 

641 



-.154 

.918 

906 




.503 

1.295 

621 



-.026 

.918 

922 




.600 

1.480 

604 



.194 

.918 

840 




.693 

1.699 

586 



.340 

.927 

750 




.781 

1.905 

578 



.507 

1.037 

650 




.856 

2.027 

567 



.591 

1.248 

626 




.949 

2.147 

566 



.657 

1.485 

609 




1.081 

2.261 

567 



.768 

1.865 

560 




1.315 

2.381 

577 



.851 

2.104 

565 




1.518 

2.457 

586 



.966 

2.226 

562 




1.760 

2.492 

589 



1.094 

2.329 

569 




2.025 

2.507 

591 



1.222 

2.396 

576 




2.541 

2.513 

591 



1.451 

2.498 

587 




3.062 

2.531 

593 



1.628 

2.560 

594 




3.490 

2.531 

593 



1.893 

2.570 

595 









2,360 

2.591 

597 









2.696 

2.560 

594 









3.123 

2.550 

593 











..J 
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APPENDIX B - Concluded 


x/d = 19.80 


x/d = 37.3 


y/d 

M 

u 

5.290 

2.494 

588 

5.043 

2.530 

592 

4.831 

2.488 

588 

4.619 

2.478 

587 

4.394 

2.473 

586 

4.195 

2.483 

587 

3.940 

2.509 

590 

3.728 

2.509 

590 

3.476 

2.499 

589 

3.172 

2.509 

590 

2.793 

2.514 

591 

2.448 

2.499 

589 

2.126 

2.488 

588 

1.893 

2.478 

587 

1.650 

2.440 

583 

1.504 

2.391 

577 

1.341 

2.324 

570 

1.240 

2.266 

564 

1.107 

2.171 

556 

-.993 

2.052 

548 

-.878 

1.870 

536 

-.821 

1.775 

532 

-.763 

1.610 

517 

-.657 

1.411 

517 

-.591 

1.303 

516 

-.490 

1.191 

522 

-.349 

1.083 

528 

-.088 

1.034 

537 

.124 

1.053 

534 

.326 

1.138 

544 

,397 

1.186 

548 

.472 

1.247 

550 

.512 

1.320 

562 

.565 

1.398 

565 

.622 

1.489 

562 

.701 

1.659 

569 

.785 

1.827 

574 

.874 

1.993 

576 

.988 

2.147 

572 

1.116 

2.249 

568 

1.271 

2.324 

571 

1.438 

2.391 

577 

1,663 

2.456 

584 

1.959 

2.488 

588 

2.387 

2.499 

589 

2.766 

2.504 

589 

3.238 

2.509 

590 

3.737 

2.494 

588 


2.236 564 
2.385 582 
2.455 590 
2,503 595 
2.501 594 
2.501 594 
2.498 594 
2.488 593 
2.488 593 
2.498 594 

2.482 592 
2.434 587 
2.357 580 
2.300 576 
2.219 571 
2.121 568 
2.052 570 
1.918 565 
1.812 560 
1.722 556 

1.618 545 
1.524 543 
1.442 540 
1.356 540 
1.340 538 
1.372 536 
1.461 535 
1.626 551 
1.760 559 
1.946 571 


-2.056 0.000 


2.507 595 
2.503 595 
| 2.488 593 


2.459 588 
2.470 590 
2.470 590 
2.470 590 
2.470 590 
2.470 590 
2.470 590 
2.470 590 
2.462 589 
2.438 586 

2.405 583 
2.361 581 
2.338 580 
2.304 578 
2.234 574 
2.155 571 


2.201 I 0.000 
1.778 


-.975 

2.036 

568 

-.856 

1.962 

566 

-.785 

1.900 

562 

-.666 

1.799 

555 

-.551 

1.717 

549 

-.335 

1.604 

544 

-.110 

1.562 

547 

.115 

1.579 

551 

.251 

1.604 

551 

.459 

1.717 

559 

.569 

1.815 

568 

.701 

1.900 

570 

.847 

2.023 

576 

.979 

2.137 

581 

1.121 

2.210 

580 

1.200 

2.257 

582 

1.381 

2.338 

585 

1.566 

2.388 

586 

1.707 

2.422 

587 

1.893 

2.449 

588 

2.113 

2.459 

588 

2,396 

2.465 

589 

2.872 

2.470 

590 

3.326 

2.476 

590 

3.834 

2.478 

590 

4.420 

2.470 

590 

4.923 

2.443 

587 

5.184 

2.394 

581 













APPENDIX C 


DATA REDUCTION DETAILS 


Primary variables measured in the test program included pitot pressures, total 
temperatures of both the hydrogen and the air jets, volumetric concentration of hydrogen, 
barometric pressure, and probe location. The hydrogen flow to the center nozzle was 
also determined by use of a calibrated sharp-edge orifice meter and redundantly mea- 
sured by a variable -area orifice meter. (The flow rates determined by the two techniques 
were in agreement within 1 percent.) 

The center of the flow field was assumed to coincide with the center of the hydrogen 
concentration (volumetric) profile. Since concentration and pitot measurements were not 
necessarily made at the same radial location, concentration values needed to determine 
velocity from the pitot-pressure measurements were taken from a plot of concentration 
as a function of probe location. The volumetric concentration values were reduced to 
hydrogen mass fraction values by assuming that the flow field consisted of a binary mix- 
ture of air and hydrogen. The mass fraction of hydrogen a is related, through the 
molecular weights, to the hydrogen volumetric concentration /3 by the following equation: 


2.016/3 

2.016/3 + 28.96(1 - $ ) 


(Cl) 


The local molecular weight is given by 


7? = 2.016/3 + 28.96(1 - /3) 


(C2) 


and the local gas constant by 



(C3) 


In equation (C3), R is the universal gas constant. (Note specific values are not assigned 
to parameters in this appendix if the values of the parameters vary with different systems 
of units. Any consistent system of units may be used in the evaluation of the equations 
presented.) The mixture specific heat Cp is expressed as 


Cp - Cp ^Q! + C p a (l - a) (C4) 

In equation (C4) , Cp ^ is the specific heat at constant pressure of hydrogen and Cp >a 
is the specific heat at constant pressure of air. The local total temperature T^ was 
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APPENDIX C - Concluded 


computed by an energy balance as given by 


T+ = 


^p,h aT tJ + c p,a( 1 _ a ) T t,a 


(C5) 


In equation (C5), T t j and a are the measured total temperatures of the hydrogen 
and airstreams, respectively. 

The local Mach number was computed from the Rayleigh pitot formula (see eq. (100) 
of ref. 18) for supersonic Mach numbers and from basic isentropic flow relations (see 
eq. (44) of ref. 18) for subsonic Mach numbers. A constant value of the specific -heat 
ratio y of 1.4 was used in all data reduction. The static pressure was assumed uniform 
and equal to local atmospheric pressure throughout the flow field. Local static tempera- 
ture was computed from the local total temperature (eq. (C5)), and the local Mach number 
was computed from isentropic flow relationships (see eq. (43) of ref. 18). 

The local velocity is then computed from 

u = M]fyRT (C6) 

In equation (C6), M is the local Mach number and T is the local static temperature. 

Equation (1) 


m x =fif 2EMdA 
Ja 

may also be evaluated to determine the total hydrogen flow rate m x at any survey 
station. In this equation p is the static pressure, and the integral is evaluated over 
the area A for which a is greater than zero. The other parameters are as previously 
defined. 
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23*114 cm (Mach 1.32 nozzle) 
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(b) Nozzle outer contour details. (c) Mach 1.32 plug nozzle contour details. 



X 


-r = 3.00 



Contour coordinates 


X 

y 

X 

y 

X 

y 
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5.930 
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5.854 
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10. 800 
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0.748 

11. 600 
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1.912 
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1 . 854 
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Note: For location of origin of coordinates relative 

to the overall nozzle arr angement, see figure 1(a) 

(d) Mach 2.50 plug nozzle contour details. 
Figure 1.- Concluded. 







Outer mixing boundaries due to outer jet 
mixing with surrounding atmospheric air 











(a) Knife edge parallel to flow 
flash exposure. 


(b) Knife edge parallel to flow 
time exposure. 


(c) Knife edge normal to flow 
flash exposure. 


(d) Knife edge normal to flow 
time exposure. 

L-71-691 

Figure 5.- Schlieren photographs for M a = 2.50 and M, = 0.91. 









(c) Air-air mixing for M a = 1.268. (d) Air-air mixing for M a = 1.302. 

Figure 7.- Center-line velocity distributions predicted by unified and Z -difference eddy viscosity 

models compared with data. 
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(c) Air-air mixing for M a = 1.268. (d) Air-air mixing for M a = 1.302. 

Figure 9.- Theoretical center-line velocity distributions as predicted by Cohen and kinematic Z-difference eddy 

viscosity models compared with data. 
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ten and kinematic Z -difference eddy viscosity models compared with data. 








Cohen viscosity model, -0.002, N 
Kinematic Z-difference viscosity mod 





(a) Velocity profiles; x/d = 0.0 to 9.58. 
Figure 11.- Hydrogen-air mixing data and correlation for 



Kinematic Z-difference 



(b) Velocity profiles; x/d = 15.4 to 63.6 
Figure 11.- Continued. 
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(c) Hydrogen mass fraction profiles; x/d = 0.0 to 9.58. 
Figure 11.- Continued. 




(d) Hydrogen mass fraction profiles; x/d = 15.4 to 63.6 
Figure 11.- Concluded. 

















(c) Hydrogen mass fraction profiles; x/d = 0.0 to 8.75. 
Figure 12. - Continued. 





data point 



Figure 12.- Concluded. 









